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SUMMARY 


The rate equation for the dissociation of air, based on simple- 
collision theory, and the numerical values of the factors that enter it 
are discussed. The rate equation and the conservation equations for a 
flowing dissociating gas are solved for the conditions (after the adjust- 
ment of vibration is taken into account) behind normal shock waves at 
Mach numbers of ID, 12, and li|-. As a result, the process of thermal 
dissociation of air is shown from its initiation to its approach to 
equilibrium. 

The results are applicable to a range of densities, and the Mach 
nixmber range covers the temperature range between a temperature below 
which little dissociation occurs and a temperature above which electronic 
excitation and ionization may have to be taken into account. The distance 
required for a moderately large fraction of the oxygen to dissociate 
varies widely with Mach nuinber and density, frcsn. a fraction of a milli- 
meter to hundreds of meters. The dissociation of nitrogen, however, is 
too slow for much dissociation to occur within a minimum distance of 
about 20 meters. It has, therefore, been established that, granted the 
theory and the numerical constants used are accurate, the dissociation 
of nitrogen outside the boundary layer cannot have any significant aero- 
dynamic effect at any altitude at Mach numbers less than l4. 


INTRODUCTION 


The shock waves that will be produced by flight at the higher hyper- 
sonic speeds will result in high air temperatures that will cause at 
least a partial dissociation of the molecules of air into atoms. This 
dissociation represents a heat sink that might conceivably alleviate the 
problem of aerodynamic heating. On the other hand, the greater thermal 
conductivity of atoms or the release of heat by the recombination of 
atoms on the surface of the vehicle might aggravate the problem. 
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Because dissociation is not an instantaneous process tut occurs at 
a finite rate^ it follows that a thorough understanding of aerodynamic 
heating at hypersonic speeds requires seme information on the rates of 
dissociation and of recombination. The preset paper is concerned with 
calculations of the temporal and spatial rates of dissociation following 
normal shoch waves and of the approach to dissociative equilibrium. 

Other available calculations with regard to the rate of dissociation 
of air are those presented by Bethe and Teller (ref. l). The present 
calculations differ from Bethe and Teller’s in a number of respects. In 
the present paper, the dissociation process is followed frem its initi- 
ation to its approach to equilibrl-um; whereas in reference 1 , only the 
rates at various temperatures and the final conditions are given. Bethe 
and Teller do not show as such the variation of the degree of dissociation 
or the variation of the state variables with time or distance behind the 
shock wave. In the present paper, furthermore, the results of the calcu- 
lations are so presented that they apply to a range of air densities j 
whereas in reference 1 , only one density was used, and it was not very 
representative of the high altitudes at which hypersonic flight must 
occur. Other differences are that a somewhat different equation is used 
herein for the rate of dissociation, an equation is derived for the rate 
of recembination, and a significantly different value is used for the 
dissociation energy of nitrogen. 

The present calculations are based on Simple-collision theory. 

Because certain numerical factors cannot be evaluated by this theory but 
must be determined by experiment, and because no experiments on oxygen 
or nitrogen have been reported, values of these factors have been used 
that were based on consideration of the available experimental results 
with other gases. Because of the lack of direct experiment and because 
rates of dissociation may depend to some extent on quantum mechanical 
effects that are not taken into account bythe simple-collision theory, 
the results of the present calculations should be considered to be uncer- 
tain within perhaps two orders of magnitude. This large degree of uncer- 
tainty is, however, in some cases not important - for example, in the 
case in which the calculations show that a significant degree of dissoci- 
ation does not occur until the air has moved 1(P meters behind the shock 
wave. 


The calculations have been made for normal shock waves at Mach num- 
bers of 10 , 12 , and l 4 , which cover most of the range between two temper- 
atvires. Below one temperature there is little dissociation, and above 
the other, electronic excitation and ionization may have to be taken into 
account. Ihe free-stream temperature for all conditions was taken as 
500° K. ■ - 


Because adjustment of the specific heats of air to vibration of the 
molecules, following a thermal jimip, is a much faster process than the 



HACA IN 363 k- 


5 


process of dissociation, the effects of this adjustment on the properties 
of the air and on the variables of state were calculated and taken into 
account as thovigh they were complete before the process of dissociation 
began. 


SYMBOLS 


a 

A,B 

D 

F° 

k 




speed of sound, cm/sec 
activation energies, ergs /mole 
constants in the conservation equations 

dissociation energy, ergs /mole 
free energy, ergs/mole 

Boltzmann's constant, I .58 X lO"^ ergs — 

molecule-°K 


3 

specific rate constant for dissociation, SoT e ^T e^ 


sec 


Kc 

% 

M 

M 




Nr 


specific rate constant for recombination, 

equilibrium constant, molecules/cm5 

equilibrium constant, atm 

molecular weight, gm/mole 

average molecular weight, gm/mole 

Loschmidt *s number, 2.687 X 10^^ molecules 

atm-cm-^ 

Mach n\miber 


Avogadro's number, 6.02 X 10^^ BP . ikicules 

mole 


cm^ 


molecule‘'^-sec 


pressure, dynes/cm^ 


la 
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P 

R 

t 

T 

u 

U 

Z 

3 

P* 

7 

7 

M- 

P 

a 

T 

( ) 


steric factor 

gas constant, 8. 31 X 10*^ 

mole-'-TC 

time, sec 
temperatiire, °K 
velocity, cm/sec 
unspecified molecule or atom 

colUslon rate, 

cm-^-sec 

energy content per molecule divided by kT 
energy content per gram divided by p/p 
ratio of specific heats 
average value of ratio of specific heats 
reduced molecular weight, gm/mole 

mass density, gm/cm^ 
molecular diameter, cm 
duration of collision, sec 
concentration, molecules /cm? 


Subscripts : 


1 free-stream station 

2 station behind shock wave where translation and rotation 

are, in equlllbrlxmi 

5 station behind shock wave where vibration is in equilibrium 

4 station behind shock wave where dissociation of oxygen is 

in equilibrium 
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5 station ‘behind shock wave where dissociation of nitrogen 

is in equllihrium 

o initial, except for Avogadro’s numher and Loschmldt’s numher 


THEORY 


When air passes through a strong shock wave, there follows a conver- 
sion of energy of directed translational motion into other manifestations 
of energy, such as randcm translation, rotation, vibration, and dissoci- 
ation. Because the temperature of a gas depends on the energy of random 
translation, when energy is transferred into this manifestation, the 
temperattnre of the gas rises. If, at a subsequent time, a portion of the 
energy of random translation is converted into other modes, the temper- 
ature falls. Because of differences in the magnitudes of the time 
required for the various processes to absorb their equilibrium amounts 
of energy, considering them to occur separately and consecutively results 
in a good approximation. Ihe variation of temperatiare frcm its free- 
stream value to its value following establishment of thermal equilibrium 
can, therefore, be represented as in the following schematic diagram, 
in which different (and nonlinear) scales are used between different 
stations : 



Station (l) represents conditions in the free stream ahead of the shock 
wave where the temperature is low enough (500° K for the present calcu- 
lations) that vibration of the molecules is not appreciably excited. 
Rotation of the molecules is, however, ccrapletely excited, and the value 
of the ratio of the specific heats 7 is 1.4. Station (2) represents 
conditions behind but very close to the shock wave where translation 
and rotation are again in equilibrium, but vibration is just beginning 
to be excited. At station (3), vibration is fully established, but 
dissociation to a first approximation is only beginning. At station (4), 
sufficient time has elapsed for the dissociation of oxygen to reach 
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equlUbriumj and at station (5)^ the dissociation of nitrogen is in 
equilibrium and the temperature has reached its equilibrium value. 


Adjustment of Translation and Rotation 

Station (2) represents conditions behind but very close to the shock 
wave. The thermal Jump across the shock wave throws the translational 
and rotational degrees of freedom temporarily out of equilibrium, but 
only a few collisions are sufficient to establish equilibrium not only 
of the translational but also of the rotational degrees of freedom 
(ref. 1). Station ( 2 ) then is taken to be' only a few mean free paths 
behind the front of the shock wave. Although the temperature at sta- 
tion (2) is high enough for the vibrational degrees of freedom to be 
excit^, sufficient time (or a sufficient number of collisions) has not 
passed for them to be appreciably excited, and the ratio of the specific 
heats remains at l.k. 


Conditions at station (2) are, therefore, calculable by the Ranklne- 
Hugoniot equations for a normal shock wave as 


>^2 1 + 0-2%a,l^ 


^ 2 
6 

%a,2^ 


%a,l^ “ I 

= %a.l^ + ^ 
7%a,l^ - 1 


> 


( 1 ) 


Adjustment of Vibration 

Adjustment of the vibrational degrees of freedom' occurs between 
stations (2) and ( 5 )* This process requires a n-umber of collisions (or 
a time or distance) that is large compared with that required for adjust- 
ment of translation and rotation, but is small compared With that required 
for establishment of dissociative equilibrium. As discussed extensively 
in reference 1 , the number of collisions required for the establishment 
of vibrational equilibrium of oxygen and nitrogen molecules is very 
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difficult to establish accurately on the basis of the current state of 
theory and experiment. The estimated values are such (lOO at high temper- 
atures), however, that the time required is much less than that required 
for dissociation, and neither the time nor the distance required is dis- 
cussed further in the present paper. The conditions at station ( 5 ) must 
be calculated, however, before calculations can be made of the rate of 
dissociation. 

The change in conditions between stations ( 2 ) and ( 5 ) is herein 
considered to occur in a single step, and conditions at station (5) are 
calculated from those at station (2) by means of the three conservation 
equations. At station ( 2 ), five degrees of freedom f are active, and 
the ratio of specific heats is 


r = ^4^ = 1.4 ( 2 ) 

Since, at this station, dp‘/^^ is zero, 

(3) 

At station ( 3 ), if there were no effects of electronic excitation, 
anharmonlcity of vibration, interaction between rotation and vibration, 
and stretching due to the centrifugal effects of rotation, the value 
of p' would be 4 . 5 * Reference 2 shows, however, that, over the range 
of touperattire Involved herein, the average value of p* for the nitro- 
gen molecule is about If . 6 and for the oxygen molecule is about if. 9 * 
Accordingly, a value of if. 7 ^or air was adopted for use hefeln. 


The equations for conservation of energy, mass, and momentvim between 
stations (2) and (3) are then, respectively. 


VI 

(M 

+ 

P2 ^2 

(1^) 

P3U3 


(5) 

PjU^S 


(6) 


These three eq^uations can be solved for p^. From equations (if) and (5)^ 

+ 0.5c2^ = ^^1^3^ 


( 7 ) 
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Frcm equations (5) and (6), 

i4-*7P3p3 + 1 »-*Tc 2^ = l4-*7c3p3 

From equations (7) and (8), 

C]^P3^ - 4.7C3P3 + 4 . 2 c 2 ^ = 0 


and 


P3 = 


^.7=3 t •7c^)^ - l6.8c^c 


2c. 


( 8 ) 


(9) 


(Only the positive sign should he used ahead of the root in order that 
the second Law of thermodynamics not he violated.) 

The velocity u^ is then given hy equation ( 5 ) and the pressure 
P 3 hy equation ( 6 ) or equation (It-). The temperature T 3 is given hy 
the general gas law 

^5 "I ‘^5^5 ■" 

in which the molecular weight M is 28.8 gm/mole. 


Rate Equations for Dissociation 

The principal purpose of the present paper is to calculate conditions 
(as fxmctions of time and distance) between stations (5) and (5). In 
order to make these calculations, it is necessary to have (l) an expres- 
sion for the rate at which dissociation occurs, (2) an expression for the 
rate at which reccmhinatlon occurs, and ( 5 ) equations that express con- 
servation of energy, mass, and momentum. The present section is devoted 
to a derivation of the rate equation for dissociation and to a discussion 
of the choice of ntimerlcal values for the quantities that appear in that 
equation. 

The simple physical picture of dissociation that is employed herein 
is that the gas is composed of diatomic molecules, that the two atoms 
of each molecule are held together hy a rather strong binding force, and 
that at station (l), ahead of the shock wave, the temperature of the gas 
is low and the kinetic energy of random motion of the molecule is small 
enough that only a very negligible portion of the molecules is dissociated. 
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At station (5); where the temperature of the gas has been raised consid- 
erably by passage through a shock wave, the energy of random translational 
motion has been increased enough that a nonnegligible degree of dissoci- 
ation is possible, because in an appreciable fraction of the collisions 
the kinetic energy involved is great enough that a collision can result 
in breaking of the binding force and in the separation of a molecule into 
two atoms. 

The present calculations of the rate of dissociation have been based 
on simple-collision theory. The only other theory that might be consid- 
ered is the activated-complex theory of absolute reaction rates of IJyring 
and others (ref. 3 )- Use of this theory, however, was not practicable 
becaiise the potential-energy stirfaces that are required apparently have 
not yet been calculated for oxygen or nitrogen. A fundamental weakness 
of the collision theory is that ccmplicated rate processes that involve 
important quantum mechanical effects cannot accurately be described or 
predicted except by quantim mechanical or wave mechanical theory. Whether 
the dissociation of oxygen and nitrogen involves Important quantum mechani- 
cal effects is not yet known. At any rate, these effects, if present, 
would not be taken into account by the single-collision theory. The 
results of calculations of rates by this theory can, therefore, not be 
depended upon to be accurate within some factor that may be as large as 
one or two orders of magnitude. 


The collision theory postulates that the rate at which molecules 
are dissociated into atoms depends on the rate at which collisions occur 
in which sufficient energy is transferred. The rate of dissociation is 
custcmarily calculated as the product of three factors; (l) the rate at 
which binary collisions occur, (2) the fraction of these collislonB in 
which sufficient energy is involved, and (5) the fraction of collisions 
involving sufficient energy that actually results in dissociation. 


Hie first of these factors, the rate' at which binary collisions 
occur,, need not be derived herein, as the deviation is straightforward 
and is given in several textbooks on the kinetic theory of gases, on 

statistical mechanics, and on gaseous reactions. The rate, in 

cm^-sec 

at which binary collisions occur between Og and U molecules, for 
example, is 


Zn. = 


_ 2 


°2 ( 1 , 2 ) 


(02)(U) 


°u' 


( 2 jtRT) 


1/2 

1 ^ 02 % , 


1/2 


( 11 ) 


The symmetry factor (l, 2 ) is equal to 1 if the two colliding mole- 
cules are of different species and is equal to 2 if they are of the same 
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species. Immediately beMnd the shock wave, 79 percent of the U. mole- 
cules are nitrogen molecules and 21 percent are oxygen molecules . As 
dissociation proceeds, the percenta-ge of molecules decreases (oxygen 
decreases faster than nitrogen) and the percentage of atcms increases . 
The diameters and the molecular weights of the atoms differ from those 
of the molecules. The procedure adopted herein is not to write average 
values of concentration, diameter, and molecular weight for each stage 
of the process but to Include in the equation a separate term for each 
constitutent of the gas. Furthermore, inasmuch as rates of collision 
of both oxygen molecules and nitrogen molecules with all species of 
particles are needed, separate equations are written for oxygen and for 
nitrogen: 


ZO2 ' + 
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The first of several difficulties or uncertainties to be encountered 
in the present calculations is now apparent; the exact numerical values 
to be assigned to the molecular diameters are not Imown. The kinetic 
theory of gases enables the effective collision diameter for the various 
transport phenomena to be calculated from experimental meas\irements ; for 
transfer of momentimi, from measurements of viscosity; for transfer of 
energy, from measurements of thermal conductivity; and for transfer of 
mass, from measurements of rate of diffusion. The results by various 
methods do not agree exactly, but usually agree within approximately 
10 percent. Kennard (ref. 4) regards the results from viscosity as prob- 
ably the most reliable data with regard to collision diameters. For 
this and other reasons to be given presently, the values of a obtained 
by this method for molecular oxygen and nitrogen have been used in the 
present investigation. The values of a taken from Chapman and Cowling 
(ref. 5 ) are 


o-Qg = 3.6 X 10“^ centimeter 
cru^ = 5.8 X 10"® centimeter 

The data listed by Chapman and Cowling also show that the following 
trends are approximately true: For a monatomic gas, a is about 75 per- 

cent as large as the a for a diatomic gas of the same molecular weight; 
and for a monatomic gas, cr is about 75 percent as large as the a for 
another monatomic gas of twice the molecular weight. On the basis of 
these trends, the molecular diameters for oxygen and nitrogen were 
chosen for use here as 


cTq = 2.0 X 10"® centimeter 
Cjj = 2.1 X 10"® centimeter 

The following reservations with regard to the choices of cr should 
be noted. In the first place, the use of a value of a calculated from 
measurements of viscosity, which is a mamentum-transfer process, can 
scarcely be wholly Justified in energy-transfer processes. In the second 
place, a may be a function of relative velocities and, consequently, 
of temperature. However, the function is unknown because the exact vari- 
ation of molecular forces with distance that should be used in calculating 
a frcm meas’urements of viscosity is not known. (The model used in ref. 5 
for calculating ff was a smooth, rigid, elastic sphere . ) One might con- 
sider, for example, using Sutherland’s model for calculating cr, but one 
could not really wholly justify such a refinement. 
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The second of the factors that make up the expression for the rate 
of dissociation is the fraction of collisions in which sufficient energy 
is involved. (Sufficient energy means energy equal to or greater than 
the dissociation energy D. ) A diatomic molecule has energy of translation, 
rotation, and vibration. Dissociation occurs when the vibrational energy 
exceeds a certain value. However, the simple-collision theory does not 
tell which of the three forms of energy (or what combination of them) is 
Important or effective in regard to transfer of energy to vibration of one 
of the two colliding molecules. In fact, it is custcmary (as is done 
herein) to Ignore the rotational and the vibrational energies of the col- 
liding molecules and to consider only their translational energy. Even 
so, another difficulty remains. Two postulates are still clearly possible; 
one, that the total kinetic energy due to the relative motion of the two 
molecules shall exceed a specified value D, and the other, that the kinetic 
energy due only to the relative motion along the line of centers at the 
moment of impact shall exceed D. The second of these postulates is the 
one that is usually used in calculations of gaseous reactions. As many 

textbooks show, the Boltzmann factor e“D/RT gives the fraction of all 
binary collisions in which D is equal to or less than the kinetic energy 
of translation due to the component of relative motion that lies along 
the direction of the line of centers at the moment of collision. This is 
the factor that is used in the present calculations. The first postulate 

leads to the factor l)e"®/^^ as the fraction of collisions in 

\RT / 

which the energy of translation is sufficient to cause dissociation. Refer- 
ence 1 used this factor with the ”l" emitted as a simplification. 

There are two reasons why the simple Boltzmann factor was 

used in the present investigation, in addition to the fact that it is 
usually used. One reason is that, when the relative veloci-ty of trans- 
lation of two molecules is broken down into two components, one along the 
line of centers at the moment of collision and the other perpendicular 
to that direction, it appears that the velocity perpendicular to the line 
of centers affects principally the angular momentum, or mutual rotation, 
of the two molecules; consequently, only the other component of velocity 
might affect the vibration of one of the molecules. 

The other reason for the use of the Boltzmann factor is now to be 
developed. Usually, or perhaps always, an additional condition is made 
that there is only a probability P and not a certainty that dissociation 
will occur in the fraction of collisions in which the energy requirement 
is met. This condition is based on the assumption that, in order for the 
translational energy- of one of the colliding molecules to be transferred 
into vibrational energy, some of the possible relative orientations of 
the two molecules are more effective than others. This probability P 
is, therefore, called the sterlc factor. There is no way of estimating 
P on the basis of simple-collision theory, so reliance must be placed in 
experiment. No experiments, however, have been reported on the thermal 
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dissociation of oxygen or nitrogen. As a matter of fact^ experiments do 

not actually give a value for P "but only for the product Po^. Refer- 
ence 6 states that experimental data on the rates of approximately six 
homogeneous^ htmolecular, second-order, gaseous reactions can he corre- 
lated by the simple-collision theory if the second of the two postulates 
about kinetic energy is accepted (reaction occurs only in those collisions 
in which the kinetic energy due to relative motion along the line of 
centers equals or exceeds D), if molecular diameters of the order of a 
few angstroms (that is, of the order of those obtained from viscosity 
measurements) are used, and if the steric factor is taken close to unity. 

Furthermore, reference 7 (PP* IT® and ITl) shows that, for several 
bimolecular reactions, the steric fantor P has a value given approxi- 
mately by the ratio of the observed rate of reaction a to the calcu- 
lated rate of reaction b when the calculations are based on the Bolt zm a n n 
factor, and where 


5 X 10^^ < a < 5 X 10 ^^ 

1 X 10 ^^ < b < 6 X 10 ^^ 

A rough approximation to the average value of P frcm, this information 
is between l /5 and 1 . For the present calculations, therefore, the 
Boltzmann fantor was used, the molecular diameters based on viscosity 
measurements were used, and the value of P was taken as l/ 2 . 

A numerical comparison of the values used in this paper with those 
of reference 1 might be of interest. Ihe present considerations have 
the rate of dissociation of oxygen proportional to P = 0 . 5 * Bethe and 
Teller have the rate proportional to PD/RT where P = 0 . 01 . 



p 

pd/rt 

°K 

(present paper) 

(ref. 1) 

4,000 

0.5 

0.15 

8,000 

.5 

.08 


The comparison shows that the rate used in the present calculations is 
several times greater than that used in reference 1. 

The numerical values of the dissociation energies that were chosen 
for the present calculations are the latest available and are probably 
the best values. The dissociation energy of oxygen has long been accepted 
as 117,200 calories per mole. A better value now appears to be 117,960 calo- 
ries per mole (ref. 8), and that value was used herein. The dissociation 
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energy of nitrogen has long been imcertain. In reference a value 
of 182^ CXX) calories per mole was used. Recent evidence (refs. 9 "to 12) 
indicates that the most probable value is 225^100 calories per mole, 
and that value was used in the pres ©at calculations. 


The rate at which dissociation occurs, in 


molecules 


cm^-sec 

the product of the rate at which binary collisions occur 
ratio of dissociations to collisions Pe“D/^T. Thus, 


^{° 2 ) - z 

dt “ ^02^® ' 


and 


-1-^ = Zig Pe 
dt ^2 




is given by 
Z and the 


w 


(15) 


If the values of Z as given by equations ( 12 ) and (15) are used, numer- 
ical values are substituted for the various constants, and arithmetical 
simplifications are made, the rate equations (1^4-) and (15) become 


- = [3.25(02) + i^.79(0) + 7.1(N2) + 

5.19(N^1.15 X 10"^(02)\(t’ e"59i4-00/T 

and 

“ "dF^ " [7.1(02) + 5.27(0) + 5.85 (N2) + 

5.7l(N^l.l5 X 10 "^^ (N2)^ g-ll55it-0/T 


Rate Equations for Recombination 

The purpose of the present section is to derive equations for the 
rate of recombination of atoms into molecules. These equations will be 
used for comparing the ratio of the rates of dissociation and recombi- 
nation with the equilibrium constant. 
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When two atonis couibine^ a stable molecule is formed, only if a portion 
of the binding energy is removed by a third body and if momentum is con- 
served. In other words, triple collisions are necessary for the recombi- 
nation of atoms into diatomic molecules. In order to calculate the rate 
at which recombination of atoms occurs, the postulate is made that recom- 
bination can occur in two ways. Two atcms can collide and if, during 
their collision, a third body collides with them, recombination may occur. 
Or, an atom and a third body can collide and if, during their collision, 
another atom collides with them, recaabination of the two atoms may occur. 

Thp method of calculation used herein is based on the idea that 
triple collisions can be considered as binary collisions between, on the 
one hand, a colliding pair of particles and, on the other hand, a 'third 
particle. An equation similar to equation (ll) can be written for rate 
of collisions. Then, multiplication by the ratio of recombinations to 
collisions (product of steric factor and Boltzmann factor) gives as the 

rate in recombination of oxygen atoms, for exiample, 

cm^-sec 


Here, (OO) represents the concentration of pairs of 0 atoms that are 
in the act of colliding with each other; and (OU), the concentration of 
0 and U particles that are in the act of colliding with each other. 

The steric factor I’qq.u ^ a measure of the efficiency of U as -a 

remover of energy and momentum from a pair of 0 atoms, and Pqu-.q 

a measure of the efficiency of U in removing energy and momentum when 
U is already in temporary collision with an 0 atom and the pair is hit 
by another 0 atom. Ohe quantities A and B are activation energies:, 
the reaction will not occur if the energies involved do not equal or 
exceed the activation energy. 

Of course, equation (l8) is useless until the concentrations of 
colliding particles (OO) and (OU) are determined. These concentrations 
will be considered as equal to the rate at which the collisions occur 
multiplied by the mean duration of a collision. Thus, 

(OO) = TqqZqq 
(OU) = 
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By use of an equation similar to equation (ll)^ the following equations 
result: 


and 


(OO) = TQQ(0)(0)aQ^(2nRT)^^^^j^j 


(OU) = 2Tqu(6)(u) 



(19) 


( 20 ) 


These. equations, are considered to give the average concentration of 
particles engaged in hlnary collisions at any given instant. 

When (OO) and (OU) as given by equations (19) and (20 ) are substi- 
tuted into equation (l 8 ), the result is 


d(02) 

dt 








OU 




^OU-O'^OU® 


■B/^T 


( 21 ) 


A similar equation can be written for the rate of recombination of 
nitrogen atoms. Subsequent to the derivation of these equations, it was 
learned that Herzfeld (ref. 15) had used' the same concept. Apparently, 
it has been customary usually to modify Herzfeld ’s ideas and to consider ■ 
collisions between the third body and each of the two particles that are 
in temporary collision with each other rather than to consider collisions 
between the third body and the complex as a whole which has been formed 
by the two colliding particles. Possibly, the reason has been that it 
was thought easier to estimate the collision diameter of an atom than 
that of two colliding atoms . This usual modification has not been fol- 
lowed herein; this is one reason that the equations used herein are 
slightly different from those given in the literature. (See for exam- 
ple, refs. 7 (p. 177 ), (p. 78 ), 15 , and 16 . ) 

Still unevaluated are the activation energies A and B, the steric 
factor P, and the collision duration t-. Because the potential-energy 
curve for oxygen, which is well established, indicates that A is zero 
for oxygen, and because no reason is Imown why they should be large enough 
to have a significant effect, A and B are assumed to be zero. On the 
basis of experiments on the recombination of bromine atoms discussed in 
reference l4 (pp. I 77 and 248), the steric factors for that reaction 
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appear to be of the order of 0.01. On the haeis of experiments on recom- 
bination discussed in reference 7 (pp. 179 sJii l80)j the steric factors 
appear to have values much closer to unity. The experiments on the recom- 
bination of iodine atons reported in reference 16 indicate average values 
for that reaction of approximately O.25. For the present study, a value 
of 0.1 was chosen for both Fqq.u and 

The mean duration of a collision can be defined as twice the quotient 
of the range of the Intermolecular forces and the mean speed of random 
translation. In accordance with the discussion in reference 1, a value 

of 3 X 10“^ centimeter has been used herein for the range of the forces. 

The average speed is shown in the textbooks to be (8RT/atn)^/^ cm/sec. 
Therefore , 


- 6 X 




1/2 


seconds 


( 22 ) 


No attempt is made herein to defend the apparent Inconsistency of using 
a much larger value for the molecular radius ct/ 2 than for the range of 
Intermolecular forces. The fact that t as given by equation (22) prob- 
ably has about the right order of magnitude (l0“^^ to 10"^^ second) can, 
however, be Indicated by a different method. Consider that two atoms 
approach each other and that their potential energy follows the curve 
shown in the sketch: 



After one vibration, the atoms will again separate, unless their excess 
energy is removed by a third body. The duration of their conjxmctlon 

then cannot exceed one period of vibration, which is about 10“^^ to 

10“^^ second. It should be emphasized, however, that the value given 
by equation (22) is an approximation and can probably be wrong by an 
order of magnitude. 
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The effective diameters of the temporary particles formed hy the 
binary collisions vere taken as 88 percent of the sum of the diameters 
of the two colliding particles in order to give diameters that were equal 
to those of oxygen and nitrogen molecules. Thus^ 

Oqq = 3*6 X 10"® centimeter 
dQOg “ 5*0 X 10"® centimeter 

Q 

Cojf = 3-T X 1-0" centimeter 
^0N2 “ 5*2 X 10"® centimeter 
cTjjOg = 5-1 10"° centimeter 

cTjjjj = 5-8 X 10"® centimeter 


O 

= 5*5 X 10“ centimeter 

When substitution is made In the rate equation (21) for t from 
equation (22), substitutions of numerical values are made for the con- 
stants, and numerical simplifications are made, then the equation for 

the rate of recombination of oxygen atoms in -becomes the 

cm-5-sec 

following equation: 

= 1.73 X 10"5V/2(0)2|l7.l(02) +’2.4(0) + 19.5(N2) + 7.^(N^ 

(25) 

Equation (24) is a similar equation, similarly derived, for the recom- 
bination rate of nitrogen atoms: 


= 1.73 X 10’^V/^(N)^jl9.8(02) + 8.6(0) + 22.6(N2) + 5 . 1 (H^ 


( 24 ) 
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Comparison of Equilibrium Constants 

A partial deck on the rate equations can be made by comparing the 
value of the ratio of the specific rate constants with the value of the 
equilibrium constant as calculated by statistical mechanics. 

A specific rate constant for the dissociation of oxygen can be 

defined by the equation 

- ^ 

A specific rate constant k^. for the recombination of oxygen can be 
defined by the equation 

^ = kr{0)2(U) 

If, for simplicity, only pure oxygen is considered, equations (16) and 
(23) Indicate, respectively, that each rate can be written as the sum 
of two rates. Thus, 


kd(°2)W - kd,l(°2f * J%,2p2)(0) 

and 


kr(of (U) = kj,^i(0)^(02) + k^^gCo)^ 

When the degree of dissociation is large, the first term on the 
right-hand side of each of the two preceding equations is small can5>ared 
with the second term. Therefore, to a good approximation, 

% = ^,2 
^r ^r,2 

At equilibrium, moreover, the rates of dissociation and recombination 
are equal, and, accordingly, 

% = Lof_ 

kr (O2) 
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But (o)^y^(02) is just the equilibrium constant Kq. The more commonly 
used equilibrium constant is related to Kc by the equation 


Kn = 


_ T 


2T5no 


Kc 


Therefore, 

yr = g ^<i,2 

^ 273no 


and, by use of equations (16) and (23), 

^ I8IT 
^ g594oo/T 

Thermodynamics shows that = e”^ The change in free-energy 

AF° can be calculated from spectroscopic information by means of statis- 
tical mechanics. An independent partial check, of sorts, on the presently 
derived rate equations is, therefore, possible by ccmparlng the value 
of Kp obtained from the calculated change in free energy and tabulated 
in reference I7 with the value of Kp obtained from the ratio of the 
specific rate constants. A comparison is made in table I. (The values 
taken from ref. I7 were adjusted to take into account the different 
value of the dissociation energy of oxygen used herein. ) 

The last column of the table shows that the ratio of the specific 
rate constants for oxygen needs to be increased by a factor of 6:1 to 8:1 
for agreement with the equilibrium constant, and for nitrogen by a fac- 
tor of 9:1- The factor is actually smaller than might be expected. In 
reference I8, for example, a factor of about 100 was needed in the ratio 
of the rate of dissociation of iodine and of bromine gas calculated by 
collision theory to the rate of recombination determined experimentally 
by E. Rabinowitch end ¥. C. Wood for the equilibrium constant so deter- 
mined to. agree with the known values of th^ equilibrium constant. In 
reference I9, a similar comparison was made in which new data on the rate 
of recombination of iodine atoms were used, and a similar factor (of 
about 100 or more) was found to be needed for agreement. (These results 
may possibly serve as additional evidence that the steric factor used 
in ref. 1 was much too small. ) 

Although the factor of 6:1 to 9:1 that is required in the ratio of 
the present rate constants may appear to be reasonably small, the fact 
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should he borne in mind that this factor does not necessarily represent 
the uncertainty in the rate constants themselves or in equations (iS), 
(ij), (23)^ and (24). Hie rate constants and the rate equations may be 
in error by much larger factors, the errors being partially canceled when 
the ratio is taken. 


Flow Equations for Dissociation Propess 

At station ( 3 ) (see sketch on p. 5 )> * 4 ie dissociation process begins, 
and at station (5) it has closely approached equilibrium. Although it 
can be shown that at equilibrium the state of a gas does not depend on 
the rate of approach to equilibrium (see ref. l), at any station x 
between stations ( 3 ) and (5) the state of the gas does depend on the 
degree of dissociation at that station. Because the gas is assumed to 
be flowing and not stationary, calculations of its state at any station x 
are to be made not only by use of the rate equations (I6), ( 1 ?)^ ( 23 )^ 
and •(24), but also by use of flow equations based on conservation of mass, 
momentum, and energy. Such equations are developed in the present section. 

The general gas law for partially dissociated air is 



where 

_ _ 28(1^2) + i 4 (n) + 52(02) + 16(0) 
(^2) + (n) + (02) + (0) 

and 


(25) 

(26) 


P = 


Nr 


28 (N2) + 14 (N) 



( 27 ) 


The equation that expresses conservation of the ratio of number of 
nitrogen atoms to number of oxygen atcms, whether the atoms are free or 
ccmbined, is 


2(N2) + (N) = §[2(02) + (0^ 


(28) 


Conservation of mass gives 


pu = p^u^ = Cg 


(29) 
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The Integrated fom of the equation of motion is 

p + pu2 = = Cj (50) 

The equation for conservation of energy is not so easily written 
down but must be derived, as is done in the following discussion. For 
example, if is the energy content of a molecule of nitrogen, 

the internal energy in the nitrogen molecules in unit volume of air is 
3jj^kT^N2) ajid the total internal -energy content of air is, in ergs/cm^. 



+ + Po 2 (° 2 ) + 


(The notation used herein differs from that in reference 2 . The quantity 
as used herein, has the same meaning as ^(^2^ in ref. 2 .J As 

ordinarily conceived, the equation for conservation of energy is written 
in terms of \anlt mass of fluid. Division of the preceding expression 
by p then gives the internal energy content per gram of air: 


^02(^2) 


28 (N2) + iHs ) + 52(02) + 16(0) 


The sum of the internal energy of a gram of fliild and its kinetic energy 
of directed motion u^/2 must be constant at all station, except for 
the energy that goes into dissociating the gas molecules, and provided 
that the external work done on the gas is included in the "internal 
energy" and the work of dissociation is not so Included. The energy 
that is absorbed in dissociating the nitrogen in a unit mass of gas can 
be seen to be 


% 2 (J» 

2P»o 

With this dissociation energy taken into account and with the aforementioned 
specification of the meaning of p, the energy balance equation is 


R'^|PN2(^a) P02(°2) + _ ^2 _ 

28 (Ng) + + 52 ( 0 ^ + 16(0) 2 


%2(N) + Doo(O) 


+ Pq2(°2)^ ^ U3^ _ 


= _2 


2 28 (N2) + 1A(N) + 32(02) + 16(0) 28(N2)^ + 52(02) 


= Rci^ 


(31) 
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This equation is vritten as though the quantities p are not functions 
of temperature. Reference 2 has a listing of the values of the various 
quantities p which shows that considering them constant gives a good 
approximation over the range of temperatures considered herein. The 
average values chosen for the present calculations are 

%2 “ 

% - 2-7 

^02 “ “-S 
h - 2-6 

Furthermore, the values given in reference 2 for the quantities p meet 
the specifications mentioned previously: they include the external work 

but do not include the energy that is used in dissociation. (Accordingly, 
the values of p do not need to be corrected for the different value 
of Djj- used herein as compared with the value used in ref. 2. ) 


Because the rate equations contain only the concentrations and the 
temperatiire as variables (besides the Independent variable t), and 
because the flow equations must be solved simultaneously with the rate 
equations, the flow equations will be reduced to only two equations 
(the two following equations) in terms of the concentrations and the 
temperature. The equations of motion and of mass conservation can be 
combined. The quantity p is eliminated from equation (50) by means 
of the gas law ( eq. ( 23 ) ) • Then, u is eliminated by means of equa- 
tion (29). The quantities M and p are eliminated by substitution 
of equations (26) and (2j)j respectively. The result is 

2 

+ W M + 2a(N2 j ^ J>(N)t\g(02) 4- Ife(d) = “3 '5") 


When u is eliminated from equation (31 ) by means of equation (29) and 
p is eliminated by equation (27), the result is 


|Pn 2 W + Pn(^) + 302(^2) + 

28(N2) + li)-(N) + 32(02) + 16(0) 2|28(R2) + 14(N) + 52(02) + l6(0^^ 

%p(N) + Do_(o) 

± ± — = Rc 4 (3: 

2p8(W2) + 1^(N) + 52 ( 02 ) + l6(oH 
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Solution of Equations 

If equations (l6) and (23) are coniblned to give the net rate of 
dissociation of oxygen. If (Ng) Is eliminated hy means of equation (2b), 

and If equations (17) and (24) are ccmblned to give the net rate of 
dissociation of nitrogen and (N2) is again e lim i n ated, then 


^ = -[29.9(02) + 18.1(0) + 1.6(N^1.15 X 1D“^(02)WV^^^°°/^ + 

t 

[92.2(02) + 40(0) - 2.4 (N^i .73 X 10-5^(0)V/^ (34) 


and 


1 ^ = 01.4(02) + 12.4(0) + 3.8(n^ 0*76(02) + 1.88(0) - 

0.5(N)]i.15 X . 

I1o 6.5(02) + 51.-8{0) - 8A(JlTjl.75 x (35) 

If, furthermore, K2 is eliminated from equations (52) and (33) 
and numerical values are substituted for the constants, then 


1.58 X 10 -"«[ 4 . 76 ( 02 ) * 2 . 88 ( 0 ) + 0 . 5 (N^I 


6.02 X 10^^C2^ 

^ = C5 


(36) 


and 

22.2(02) + 11.2(0) + 0.4 (h) ^ 2.19 X lO^^Cg^ 

68.67 [2(02) t (og |68.67 [2(02) t (0^1" 

36,670 (n) + 29,700(0) ^ 

68.670(02) + (oj] 


(37) 
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For the actual process of coeluting the rates of dissociation, some 
simplifications were made. First, it vas known that the dissociation 
of oxygen would he shown to he essentially cranpleted "before the dissoci- 
ation of the nitrogen had hardly begun. Therefore, the ccmputational 
procedure was simplified hy computing separately the dissociation of the 
two gases. The dissociation of oxygen was carried ttirough almost to 
completion. Then, with the final conditions for that process as the 
initial conditions for the dissociation of nitrogen, the conputation of 
the dissociation process of nitrogen was carried out. 

Second, it was known that, until the dissociation was close to 
equlllhrlvun, recomhination would have an entirely negligible effect on 
the net rate of dissociation. Therefore, the computations were greatly 
simplified hy not using the last tem in equations ( 34 ) and ( 35 )* 

Furthermore, it was noticed that, because recanhination was neglected, 
the variables could be so modified that the equations would be independent 
of the initial density of the air. The rate equations, equations ( 34 ) 
and. (35 "Without their last terms were made Independent of the initial 
density by dividing them by (02)0^^ which is proportional to the square 

of the initial density; thus. 


r 

- - 29 9 + 18 1 <°) 

I (°2)o 


w. 


3.. 1,5 X 10-1^ ll/ 2 ^- 5 * 0 O/T 

( 02 ) o 


( 36 ) 


(N) 

i°2l 




21.4 ^ + 12.4 ^ + 5.8 

(°2)„ (°a)o (°s)o 




1.88 -M. . 0.5 iSl 
(02)0 W , 


1 . 15 X io-iarV 2 e-ll 554 o/T ( 39 ) 


For the present computations, (n) in equation ( 38 ) and in equation ( 39 ) 

were taken equal to zero. 
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The-flov equations (eqs. (36) and (37)), in which both the concen- 
trations Ld the constanis^ Cg and C5 (but not are proportional 

to the density, were made independent of the initial density in a similar 
manner; thus. 


1.38 X 10' 


-16 




T + 


6.02 X 10 


23 ^2 


68.67 


(°2), 



*M- 

2 

L.^ 

lo (°2)o 


(°2)o 


(ito) 


22.2 ISaL ^ 11.2 ^ OA ^ 

(O2), (°2)o (°2)= 


2.19 X 10^^ 


«.67|2 


T + 




68.67 


(=>2)0 M 


56,670 


(W) 

(° 2 )< 


+ 29,700 


(0) 

W) 


o _ 


68.67 


(Q2) . 


(02)0 (° 2 )o 


= ci^ 


(in) 


m the actual calculations, the initial conditions were taken to 
be those at station (3)- Thus, 

(02)0 ^ (^2)5 

Because of the aforanentioned modlficj-tlons to the equations, 
the results of the solution of the equations for a given value of Tq 
arP annlicable for all initial densities and can be shown by a si^le 
Tlot (?=- To) for all cleaaltl.a, with the follovln* reeehvation 

being kept in mind. The recombination terms that were omitted from 
enuatlons (38) and (39) are proportional to a higher power of “ 

centration than are the dissociation terms. Consequently, at high en- 
sl?fas! het^eeh the t»o r^as oceurs at a^Uer de^ea of 

dissociation than at low densities. Therefore, at high densities. 
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dissociation does not continue until the concentration of O2 (or of 112) 

approaches zero. This reservation, that the results as calculated can he 
expected to become inaccurate at a value of the time parameter (^2)0^ 

that depends on the initial density, is discussed further in "Results 
and Discussion." 

The computations consisted of a simultaneous solution of equations (38) 
(Ij-O), and (l 4 -l) and a simultaneous solution of equations ( 39 )j (^O), ^d 
( 4 l). The computations were made on a modified Reeves Electronic Analog 
Computer. The uncertainty in the results due to analog-ccmputer error is 
believed to be less than 5 percent. 

In order to calculate the Mach number at station ( 4 ), where the 
dissociation of oxygen is complete, the speed of sound was taken as 



where M is given by equation (26) and 7, the average specific-heat 
ratio, is 


7 = 



if the rate of change of degree of dissociation with density can be neg- 
lected, as it can be here. If the dissociation of oxygen is complete, 
then, by equation (28), 


Ship's) + 

(% - i)(%) + (Po - 4(0) 

§ PN2 + PO 


Therefore, 


7 = 1;55 


a = 2,170T^^^ cm/sec 


and 
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RESCJLOB AND DISCUSSION 

Adjustment of Translation, Rotation, and Vibration 


The results of the ccamputations of the changes in the state of the 
gas due to the adjustment of translation, rotation, and vibration behind 
the shock wave are shown in table II for three values of the free-stream 
Mach number. Conditions behind the shock wave after adjustment of mole- 
cular translation and rotation were calculated by use of equations (l)> 
Conditions after adjustment of molecul^ vibration were calculated by 
use of equations (5) ^ ^ (9)> sn^ (lO)* large drop in temperature 

due to the transfer of energy into vibration is to be noted. 

The concentration of osqygen molecules after the ccmpletion 
adjustments (02)3^ which is taken as the initial concentration 

beginning of dissociation (02 )q> calculated for a given 

stream density by the equation 

{02), " 4.38 X 10 ^^ -2. -i. p- 
which was obtained from equations (27) ani (28). 


of these 
at the 

free- 


Adjustment of Dissociation 

The results of the computations of the dissociation process are shown 
in table II and figures 1 and 2. Figure 1 shows the course of the disso- 
ciation process of oxygen for three values of the temperature at sta- 
tion (3) (three values of the free-stream Mach number). Ihe fact that the 
time rate of dissociation depends greatly on the tenperature is clear 
frcan an examination of the change in the magnitude of the horizontal 
scales in figure 1 as the initial tanperature is changed. The cooling 
of the gas due to dissociation is shown by the curves of T. The fact 
that the process takes place at nearly constant pressure is shown by 
the curves of pressure ratio. The increase in the density of the gas, 
due to cooling as dissociation proceeds, is shown by the ciirves of den- 
sity ratio. The density ratio is less, however, than the inverse of 
the temperature ratio, because the molecular-weight ratio decreases as 
the degree of dissociation increases. 

In figure l(a), especially, the degr&e of dissociation at 

6.4 X 10^^ mo lecule -sec jg greater than might appear at first glance 
cm5 

because of the density Increase. The value of (02)/(^2 )q per- 

cent, but the degree of dissociation of oxygen is not 100 - J2, or 
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28 percent; but^ because the density has increased by the ratio 1.26; the 
degree of dissociation of oxygen is LOO - 57 j or U-3 percent. 

Figure 2 shows the dissociation process of nitrogen for a free-stream 
Mach number of 14. (The conditions at t = 0 in fig. 2 were obtained by 
extrapolation of the cuirves of fig. l(c) to complete dissociation of 
oxygen. ) The much lower speed of the dissociation process of nitrogen 
compared with that of oxygen is obvious from a comparison of the horizon- 
tal scales of figures l(b) and 2 since both processes begin at approxi- 
mately the same temperature. 

The choice of variable for the horizontal coordinate of figures 1 
and 2 implies that the results shown in these figures apply at all air 
densities, as they do with the following exceptions: In obtaining the 
results shown, the recombination process was neglected. Because the rate 
of recombination (the "rate" referred to is defined subsequently in this 
paragraph) is proportional to the density and the rate of dissociation is 
independent of the density, recombination at high densities cannot be 
neglected over the entire extent of the ciirves. The upper limit of the 
density for which the results as shown are applicable can be determined 
by comparing the rates of dissociation and recombination as given by the 
equations developed in the present paper for conditions at the right-hand 
end of the curves in the figures, and by determining, the density at which 
the rate of recombination becomes as large as, for example, 10 percent of 
the rate of dissociation. In the ccmparison shown in table III, the rate 
of recombination of oxygen was divided by 6.8 for a tanperature of 6 , 560 ° K, 
8.1 for a temperature of 4,630° -K, and 8.9 for a temperature of 3,450° K; 
and the rate of recombination of nitrogen was divided by 9 in. order to make 
the ratio of the specific-rate constants agree with the equilibrium con- 
stant. (See the section entitled "Ccxnparison of Equilibrium Constants.") 
The rates given for oxygen in table III refer to the absolute value of the 
quantity 


and for nitrogen the rates refer to the absolute value of the quantity 



(K) 

VWo 



The limiting value of (Opjo ®° make the rate of recom- 

bination about one-tenth the rate of dissociation at the right-hand end 
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of the curves. The limiting density (or p^) corresponds to the 

limiting concentration ( 09 )^^ ajid the limiting density p^ is the free- 

stream density that results in the density Pq at station (5)> The 
limiting altitude, the density altitude corresponding to the limiting p^_, 
was taken from table II of reference 20. 

Another check on the upper limit of the density for which the full 
extent, of. -the curves is applicable can be obtained frcxa table 5 of refer- 
ence 17 which gives the equilibrium composition of air as a function of ^ 
pressure. The table is a rather rough approximation for the present 
application because it was calculated for equilibrium of both the oxygen 
and the nitrogen dissociations, whereas equilibrium was not reached in 
the present calculations nor, when the dissociation of oxygen was calcu- 
lated, was the dissociation of nitrogen included. Nevertheless, use of 
table 5 of reference I 7 gives estimated altitudes that are qualitatively 
in agreouent with those shown in table III of the present paper. 

The limiting altitudes shown in table III apply, as mentioned, at 
the ends of the curves (except for the bottom row of the table). The 
curves of figure 1 at smaller (less than full-rscale) values on the hori- 
zontal axis would be applicable at lower altitudes (higher densities) 
because the slope of the. (02)J(02jQ curve increases as the origin is 

approached, and the increased rate of dissociation may compensate for 
the Increased rate of recombination that is due to the increased density. 

The tentative tables for the upper 'atmosphere (ref. 20) show a change 
in the composition of air at high altitudes due to the dissociation of 
oxygen by solar radiation. The tables show complete dissociation of oxy- 
gen at altitudes above 100 kilometers in the day and above 120 kilometers 
at night, and partial dissociation down to 80 kilometer's in the day and 
105 kilometers at night. The present results should, therefore, not he 
expected to apply at altitudes high enough that a significantly large 
portion of the oxygen is originally in a dissociated state. 

Figure 1 can be used to determine the distance behind the shock_ 
wave at which a moderate or a large portion of the o^gen is dissociated. 

For example, at a Mach number of 10 with ( 02 )^^ = 6.4 X 10^^, the time t 

is 0.13 second for an altitude of 70,000 meters. For a local velocity 
h / 

of about 4 X 10 cm/sec, the corresponding distance is about 5^000 centi- 
meters. On the other hand, at a Mach number of l4 with ( 02 )Qt = 1 X 10^, 

r 

the value of t . is 5 x 10 second for an altitude of 55 j 000 meters. For 

a velocity of about 5 X 10 ^ cm/sec, the corresponding distance is 0 .I 5 centi- 
meter. These distances are, of course, in the undisturbed main stream 
behind a normal shock wave. ■ 
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Figure 2 shows that the dissociation of nitrogen is veiy slew (com- 
pared with that of oxygen) ^ and table IH shows that appreciable disso- 
ciation of nitrogen occurs only at very small densities. For small den- 
sities, however, the time for appreciable dissociation to occur is large 
(compared with the length of a vehicle divided by the velocity of the air 
behind the shock wave). These facts fairly well establish (within the 
assumptions used in the present calculations) that the dissociation of 
nitrogen for flight at Mach numbers less than 14 at any altitude cannot ■ 
proceed far enough to cause any significant effect on the properties of 
air and, therefore, will not produce any significant aerodynamic effect 
within the length of any conceivable vehicle.^ (Referring to the entire 
length of a vehicle is being conservative, as exjjanslons around corners 
would reduce the air temperature, the rate of dissociation, and the equi- 
librium degree of dissociation. Use of 500° K probably is also conserva- 
tive, as ref. 21 shows measured values considerably below that.) This 

result is substantiated by the following examples: For (Ogjo^ = 6.4 X 10^^ 

in figure 2, at the limiting altitude of 90,000 meters, (02)0 “ ^ 10^ 

and, therefore, t = I80 seconds. Even for (02)0^ “ 0.064 X 10^^, 

t = 1.8 seconds. Inasmuch as u « 5 ^ ich cm/sec, the corresponding 

distance behind the shock wave is 9 X 1(A centimeters-. At higher alti- 
tudes (higher than the limiting altitude) where (^2)0 smaller, t 

is greater. At lower altitudes, the results at ^Og^ t = 6.4 X 10^^ do 

not apply, but those at smaller values of (02)Qt may apply. Consider, 

for example, (02)0’^ “ 0.2 X 10^^. (Below this value the degree of disso- 
ciation is certainly too small to have any important aerodynamic signifi- 
cance . ) The results will be applied at the lowest altitude for which the 
recombination rate is not too large. In the bottom row of table IH for 

nitrogen, the limiting value of ( 02)0 determined to be 4.4 X 10^^. 

Therefore, t = 0.04 second and tu = 2 X 10^ centimeters. To recapitu- 
late, at high and low altitudes the dissociation of nitrogen is too slow 
at Mach numbers less than l4 for a sigsaificant degree of dissociation to 
occiar within the length of a vehicle. 

One other condition for the dissociation of nitrogen should be consid- 
ered. For flights at altitudes where the oxygen is completely dissoci- 
ated by the sun's radiation (100,000 meters and above in the daytime), 
the dissociation of nitrogen will not be preceded by the thermal dissoci- 
ation of oxygen (which causes a large drop in toi^erature) and, therefore, 
will occur at a higher temperature and at a greater rate. A simple 


4t should be stressed that this result does not necessarily apply 
for the slower flow in the boimdary layer. 
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approximate calculation is made here to determine whether the rate is 
sufficiently increased that a significant degree of dissociation can 
occur within the length of a vehicle. The only change in the rate will 
he that due to increased temperature, because in both cases complete 
dissociation of oxygen occurs before the dissociation of nitrogen begins. 
The rate will, therefore, be increased by the product of the ratio of 
the square roots of the temperattires and the ratio of the Boltzmann fac- 
, tors. In one case, the initial temperature is T^ (equal to 8 , 857 ° K), 
and, in the other, is (equal to 6,614-0° K). The initial rate is, 

therefore, Increased by a factor of 114-2. The smallest value of (*^ 2 ) 0 ^ 

that is consistent with a reasonable degree of dissociation is about 
0.2 X 10‘‘'° in figure 2. With the new rate, this value becomes approxi- 
mately 1.4 X 10 ^ 5 . At the lowest altitude (100,000 meters), the time is 
smallest ( 0.16 second). If u » 5 10^ .cm/sec, this time corresponds 

to a distance of about 8 X 10^ centimeters. 


CONCUJDING REMARKS 


The present calculations of the rate of dissociation of air have 
been based on available theory and available experimental results with 
other gases. The results of the calculations represent an extension of 
Bethe and Teller's work in that the present results follow the dis- 
sociation process frcan its Initiation to the approach to equilibrium 
and in that they apply to a range of densities . The results show that 
the distance behind a normal shock wave required for a moderately large 
fraction of the oxygen to dissociate varies widely with Mach number'and 
density, from approximately hundreds of meters to a fraction of a milli- 
meter. ,The distance reqiiired for a significant portion of the nitrogen 
to dissociate (outside the boundary layer) is shown to be greater than 
the length of a vehicle at any density and at any Mach number less than 
14. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 5, 1956. 
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OABLE I.- COMPARISON OF EQUILIBRIUM CONSTANTS 


Temperature, 
°K • 

% 

(from ratio of rate 
constants ) , atm 

(.c XT, 

corrected), atm 

Ratio of values 
of Kp 

Oxygen 

4,000 

0.2576 

2.182 

0.118 

4,500 

1.506 

12.21 

.123 

5,000 

6.270 

48.42 

.129 

5,500 

20.32 

149.8 

.156 

6,000 

54.49 

384.8 

.142 

7,000 

261.5 

1,706 

.155 

8,000 

863.2 

5,258 

.164 

Nitrogen 

4,000 

3.285 X lO""^ 

5.087 X lO"^ 

0.106 

4,500 

8.609 X 10"^ 

7.956 X 10"5 

.109 

5,000 

1.187 X 10"^ 

1.077 X 10"^ 

.110 

5,500 

1.025 X 10“5 

9.198 X 10"5 

.111 

6,000 

6.229 X 10 "^ 

5.551 X 10 

.112 

7,000 

1.080 X 10"^ 

9.629 X 10"^ 

.112 

8,000 

9.540 X 10 “^ 

8.503 

.110 




■ti£9^ VDW 









































